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ABSTRACT 

Extremely low mass (ELM) white dwarfs (WDs) with masses < 0.25 Mq are rare objects that 
result from compact binary evolution. Here, we present a targeted spectroscopic survey of ELM WD 
candidates selected by color. The survey is 71% complete and has uncovered 18 new ELM WDs. Of 
the 7 ELM WDs with follow-up observations, 6 are short-period binaries and 4 have merger times 
less than 5 Gyr. The most intriguing object, J1741-I-6526, likely has either a pulsar companion or 
a massive WD companion making the system a possible supernova Type la or .la progenitor. The 
overall ELM Survey has now identified 19 double degenerate binaries with <10 Gyr merger times. 
The significant absence of short orbital period ELM WDs at cool temperatures suggests that common 
envelope evolution creates ELM WDs directly in short period systems. At least one-third of the 
merging systems are halo objects, thus ELM WD binaries continue to form and merge in both the 
disk and the halo. 

Subject headings: binaries: close — Galaxy: stellar content — Stars: individ- 
ual: SDSS J011210.25-hl83503.7, SDSS J015213. 774-074913.9, SDSS 
J144342.74-M50938.6, SDSS J151826.68+065813.2, SDSS J174140.49+652638.7, 
SDSS J184037.78-h642312.3 — Stars: neutron — white dwarfs 



1. INTRODUCTION 

Extremely low-mass (ELM) WDs with masses <0.25 
Mq are created when the progenitor stars lose so much 
mass during their evolution that they never reach helium 
burning. Binary evolution is considered the m ost likely 
origin for low-mass WDs (e.g. lMarsh et al.ll995| ). making 
ELM WDs the signposts for the type of compact systems 
that are stro ng gravitational wa ve sources. In Paper I 
of this series (jBrown et al.ll2010i) we presented the first 
complete, well-defined sample of ELM WD s fortuitously 



targeted by the Hypervelocity Star Survey (iBrown et al 
20051 |200d |2009( ). In paper II of this series (Kili c et al 



2011al ) we characterized other ELM WDs identified in 
the Sloan Digital Sky Survey (SDSS) Data Release 4 
(|Eisenstein et al.ll2006[ ). In both programs ELM WDs 
exist in < 1 day orbital period binary systems, with an 
estimated merger rate comparable to the rate of under- 
luminous supernovae (|Brown et al.ll20l"Tal ). 

Here we present the results of the first targeted sur- 
vey for new ELM WDs. Our approach is to select high- 
probability ELM WD candidates from a well-defined re- 
gion of color and apparent magnitude, and then obtain 
spectroscopy for the objects from that selection region. 
The spectroscopic survey is presently 71% complete and 
contains 21 ELM WDs defined by 5.0 < \ogg < 7.0 
dex {g in cm s~^) . Three of the ELM WDs were pre- 
viously identified (jKilic et al.l 120091: IBrown et all I2010D 
and 18 are new discoveries. We have obtained follow-up 
spectroscopy for 7 of the new ELM WDs, and calculate 
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orbital solutions for the 6 with significant velocity vari- 
ability. Four of these new ELM WD systems have grav- 
itational wave merger times less than 5 Gyr. The most 
interesting object, J1741-I-6526, has a minimum compan- 
ion mass of 1.1 Mq. Thus J1741-I-6526 is most likely a 
pulsar binary, or, if the orbit is e dge-on, possibly a su - 
pernova Type la or .la progenitor (jBildsten et al.|[2007l) . 
It will begin mass transfer in <170 Myr. 

Our targeted ELM survey will yield a clean, non- 
kinematically selected sample of WDs. Once completed, 
we can use our sample to constrain the space density, 
period distribution, and merger rate of ELM WDs in 
double degenerate systems. SDSS, on the other hand, 
has not found large numbers of ELM WDs because they 
have not targeted them; existing SDSS WD spectroscopy 
comes from different target selection programs observed 
with different completenesses (Eiscnstcin ct al.. ,2006.) . 
In a stellar evolution context, our survey complements 
studi es of WD binaries wi t h main sequence compan- 
ions (jZorotovic et al.l 120111 : INebot Gomez-Moran et ahl 
120111 ). Because of ELM WDs' low surface gravities 
logg ^ 6 and small ~1 Rq orbital separations, a grow- 
ing number of systems exhibit some combination of 
tidal distortion, relativistic beaming , reflection effects , 
and eclipses (e.g. IBrown et al.ll201l"U iKilic et al.ll2011cl: 
Parsons et 3111201 lU Pvrza s et al.ll2011l: ISteinfadt et al.l 
2010; Venn es et al.l l2011). We expect that on-going pho- 
tometric follow-up will provide improved constraints on 
the nature and orbital inclination of our new ELM WD 
binary systems. 

We organize this paper as follows. In Section 2 we 
discuss our new survey design, observations, and data 
analysis. In Section 3 we present the orbital solutions 
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Fig . 1 . — Color-color diagram showing our ELM Survey selection 
region (blue bounded region) in the context of SDSS star counts 
(black dots). The ELM Survey selection targets a region of color 
space that contains 5 < logg < 7 dex hydrogen atmosphere WDs 
(magenta solid lines) and avoids known concentrations of A-type 
stars and normal WDs. Symbols indicate the spectroscopic iden- 
tification of ELM Survey targets: red stars are the DA WDs, red 
stars with yellow centers are the ELM WDs, black x's are nor- 
mal A-type stars, green circles are DZ WDs, cyan diamonds are 
quasars, and magenta triangles are other miscellaneous objects. 

for six new ELM WD binaries. In Section 4 we discuss 
the overall properties of the ELM WD sample and high- 
light correlations between temperature, orbital period, 
and secondary mass. We conclude in Section 5. 

2. DATA AND ANALYSIS 
2.L ELM Survey Design 

The ELM Survey is a spectroscopic survey of 15 < 
go < 20 low mass WD candidates selected by color. We 
use de- reddened, uber-calibrated point spread function 
magn itudes from SDSS Data Release 7 ()Abazaiian et al.l 
I2009D . Our color selection strategy, illustrated in Figure 
[U is constructed as follows. 

First, we target objects with colors consistent with the 
effective temperat ures of luminous ELM WDs. Updated 
IPanei et al.l (j2007D evolutionary tracks for He-core WDs 
indicate that 0.17 Mq WDs spend ~1 Gyr with lumi- 
nosities of 8 < Mp < 9 mag at temperatures ~10,000 K. 
iSerenelh et all (|2001f ) tracks give similar results. Thus we 
target —0.25 < (5 — ^)0 < 0.1. Second, we target objects 
with {u — g)o colors consistent with the 5 < logg < 7 sur- 
face gravities of ELM WDs. To make this color selection, 
we use a polynomial fit to the synthet ic photometry of 
DA WD hydrogen atmosphere models ()Koesteiil2008l) . 

{u - g)o < 0.83 - 1.074(g - r)o - 1.4939(g - r)l 
-t-0.8156(g - r)l + 33.42316(5 - r)^ 
-1-280.88439(3 - r)l - 492.62139(.g - r)^ 
-1993.9254(3 -r)^. (1) 

Finally, we restrict our color-selection to the most prob- 
able low mass WD candidates. We restrict our selection 
to those objects bluer in {u — 3)0 than the observed pop- 
ulation of A-type stars (our zeropoint in Equation 1) 



and redder in (u — g)o than the observed population of 
normal DA WDs, {u - g)o > 0.542(5 - r)o + 0.685. We 
exclude quasars based on their non-stellar color {r — i)o < 
1.8 — O.I501 a limit that becomes more restrictive at faint 
magnitudes where we expect greater contamination from 
increased photometric errors and from increased quasar 
number counts. Put together, this color selection strat- 
egy maximizes the contrast of ELM WDs with respect 
to foreground and background populations. 

Our target selection identifies 505 ELM WD candidates 
with 15 < 5 < 20 over ~10,000 deg^ of the Sloan Data 
Release 7 imaging footprint. Spectra for 116 of these can- 
didates already exist: 31 of the cand idates were observe d 
by the Hypervelocity Star survey ([Brown et al.l 120091) . 
and 85 were observed by SDSS. Three of the objects with 
existing spect ra are previously identified merging ELM 
WD binaries ()Kilic et al.ll2009i: IBrown et al.ll201CI|) . thus 
our initial expectation is that we will find about a dozen 
new merging ELM WD binaries in the full survey. There 
remain 389 candidates to observe. 

2.2. Spectroscopic Observations 

We obtained spectra for 245 of the ELM WD candi- 
dates in observing runs starting in 2010 September and 
ending in 2011 June. We observed 164 objects with 
g > 17 mag at the 6.5n i MMT telescope usi ng the Blue 
Channel spectrograph ([Schmidt et al.l1l989[ ). We oper- 
ated the Blue Channel spectrograph with the 832 line 
mm^^ grating in second order, providing wavelength cov- 
erage 3650 A to 4500 A and a spectral resolution of 1.0 - 
1.2 A, depending on whether a 1" or 1.25" slit was used. 
At 5 = 19 mag we used a 400 sec exposure time to obtain 
a signal-to- noise (S/N) of 7 per pixel in the continuum 
and a ~10 km s~^ velocity error. All objects were ob- 
served at the parallactic angle, and a comparison lamp 
exposure was obtained with every observation. 

We observed 81 objects with g < 17 mag in queue 
scheduled time at t he 1.5m FLWO telesc ope using the 
FAST spectrograph (jFabricant et al.lll998l ). We operated 
FAST with the 600 line mm~^ grating and a 2" slit, 
providing wavelength coverage 3500 A to 5500 A and a 
spectral resolution of 2.3 A. At 5 = 16 mag we used a 900 
sec exposure time to obtain a S/N of 10 per pixel in the 
continuum. A comparison lamp exposure was obtained 
with every observation. 

We process the spectra using IRAEQ in the standard 
way. We flux-calibrate using blue spectrophotometric 
standards (|Massev et al.|[l988l ). and we measure radial 
velocities using the cross-correlation package RVSAO 
()Kurtz fc Minkl ll998D . During the course of the ELM 
Survey we obtained repeat observations for 7 of the newly 
identified ELM WDs. 

2.3. Spectroscopic Identifications 

Of the 361 survey targets for which we have spec- 
troscopy: 285 (78.9%) are normal A-type stars with 
log 5 < 5, likely blue horizontal branch stars or blue 
stragglers in the halo; 23 (6.4%) are quasars; 8 (2.2%) 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 2. — Model fits (smooth red lines) overplotted on the com- 
posite observed spectra (black lines) for the six newly identified 
ELM WDs with significant velocity variability. We use the spec- 
tral continua to provide improved TcS constraints. 

are DZ WDs that show strong Can H and K hnes; 5 
objects labeled "other" in Figure [T] are an emission line 
galaxy and four featureless continuum objects; and, most 
relevant to this paper, 40 objects (11.1%) are probable 
DA WDs with logg > 5, of which 21 (5.8% of the survey) 
are probable ELM WDs with 5.0 < logg < 7.0 dex. 

2.4. Stellar Atmosphere Parameters 

We perform stellar atmosphere model fits us- 
ing an upgraded version of the code described by 
lAllende Prieto et all ()2006[ ) and synthetic DA WD spec- 
tra kindly provided by D. Koester. The grid of WD 
model atmospheres covers effective temperatures from 
6000 K to 30,000 K in steps of 500 K to 2000 K, and 
surface gravities from \ogg = 5.0 to 9.0 in steps of 0.25 
dex. The model atmospheres are calculated assuming 
local thermodynamic equilibriu m and in clude both con- 
vective and radiative transport (|Koesterir2008i ). 

We fit the full flux-calibrated spectra as well as the 
continuum-corrected Balmer line profiles. The spectral 
continuum provides improved constraints on effective 
temperature but exposes our fits to possible flux cali- 
bration problems - a concern for some of the objects 
observed in the 1.5m FAST queue, but generally not an 
issue for the g > 17 mag objects observed at the MMT. 
When we compare best-fit solutions, the flux-calibrated 
and continuum-corrected parameters differ on average by 
540 ± 700 K in Tcff and 0.06 ± 0.08 dex in log g. We take 
these differences as our systematic error. We consider 
the flux-calibrated parameters more robust, because syn- 
thetic photometry of the flux-calibrated model fits pro- 
vide better agreement with SDSS photometry in all five 
filters. 

Table [1] presents the Tcs and log 5 values for the 40 DA 
WDs identified in the ELM Survey. For the 7 newly iden- 
tified ELM WDs with multiple observations we use the 
flux-calibrated, summed spectra to derive T^s and \ogg, 
and we use the scatter of the flts to the individual spec- 
tra to derive errors. The remaining DA WDs typically 
have single-epoch, S/N~8 per pixel spectra and increased 
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Fig. 3. — Surface gravity vs. effective temperature of the ob- 
served DA WDs (solid red stars) and ELM WDs with <0.25 
Mq (stars with yellow centers) found in this survey, compared 
with predicted tracks for He WDs with 0.16-0.45 Mq (blue lines, 
Panel ct al. 2007) and CO WD s with 0-5-0. 7 Mq [green lines, 
iBergeron et a l. 1995; Holber g &: Bergeronll2003) . 

statistical uncertainties in their atmospheric parameters. 
Two of the more massive DA WDs, J084325.09-K371551.7 
and J095353.66 +410927.4, have SDSS stellar atmosphere 
measurements (|Eisenstein et al.l 120061 ) that differ from 
our measurements at the 1- to 2-a level. Figure [5] vi- 
sually compares our best-fit atmosphere models with the 
summed spectra of the six ELM WDs with radial veloc- 
ity variability. We attribute imperfect continuua fits to 
imperfect flux calibration, notably for J1741-I-6528 and 
J1840-I-6423 which were observed at high airmass. 

Figure [3] plots the T^s and \ogg of ah 40 DA 
WDs i n our targeted s urvey in re lation to the im - 
proved iPanei et al.l (|2007D tracks (see lKilic et"alll2010bf ) 
for He-core WDs wit h masses 0.16-0.45 Mq and the 
IBergeron et al.l (|1995f FI tracks for normal CO-core WDs 
with masses 0.5-0.7 Mq. The gap between the 0.17 
Mq and 0.18 Mq He-core WD tracks is hnked to the 
threshold for th ermonuclear flashes in the hydrogen shell 
burning phase (|Panei et al.ll2007l ). The presence of ob- 
jects with logg < 6 in the gap between tracks makes 
precise WD mass and luminosity esti mates diflicult (see 
iKilic et al.|[2011aHVennes et al.ll201lD . More reliable es- 
timates are possible for the log g > 6 WDs. Mass, lumi- 
nosity, and heliocentric distance estimates are presented 
in Table [H 

2.5. Radial Velocities 

We maximize our sensitivity to velocity variability with 
the following approach. We first cross-correlate the indi- 
vidual spectra with a high signal-to-noise WD template. 
We then shift the individual spectra to the rest frame, 
and sum them together to create a template for each 
object. Finally, we cross-correlate the individual spec- 
tra with the appropriate template to obtain the final ve- 
locities for each object. The average precision of our 
measurements is ±12 km s~^. We verify our velocities 

^ http:/ /www. astro. umontreal.ca/~bergeron/CoolingModels/ 
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TABLE 1 

DA White Dwarf Physical Parameters 



Object 






90 
(mag 


') 






(mag) 


Q 




{q — r)n 
fmatj) 




(K) 




log 


g 






Ma 
y 




(kpc) 


Mass 
(Mq) 


JUUiDZZ, 


no nn /I QOQ 
Uy— UU4v5Zv5. 


,4 


19, 


,704 + 


0, 


,022 


0, 


DUZ 


1 n 
iL U. 


,078 


-0, 


,239 + 0, 


,040 


11190 + 260 


7, 


A n 
4U 


it 


0, 


,08 


10 


.90 + 0, 


,14 


U. 


.oo 


U. 


,oO 


JUiiziU. 


,zoH-ioooUo. 


,7 


17, 


,111 ± 


0, 


,015 


0, 


oi n 

yiy 


1 n 
iL U. 


,025 


-0, 


170 + 0, 


,017 


9690 ± 150 


5, 


Ov5 


1 

it 


0, 


,06 


8 


.00 + 0, 


,20 


U. 


DO 


n 
U. 


1 fi 
.10 


JUli ( ZD. 


y1 0_l_0fi1 Q/1 Q 
,4yH-Z0iv54o. 


,2 


19, 


,300 + 


0, 


,018 


0, 


o 1 y 


it u. 


,055 


-0, 


,238 + 0, 


,027 


10710 ± 190 


7, 


.Do 


it 


0, 


,08 


11 


.47 ± 0, 


,15 


n 


QV 
.O 1 


n 
U. 


A fi 
,40 




,v5 ( -|~4DiyZU. 


1 


15, 


,835 ± 


0, 


,012 


0, 


. / 1 O 


it u. 


,022 


-0, 


,041 + 0, 


,014 


9050 ± 590 


7, 


1 fi 

. io 


it 


0, 


12 


11 


.30 + 0, 


,33 


n 


nfi 

.Uo 


n 


oo 

,zy 


Tm f^OI Q 


7'7'_l_n7'/101 


,9 


18, 


,015 ± 


0, 


,014 


0, 


fil O 

oiy 


it u. 


,034 


-0, 


,249 + 0, 


,021 


10840 ± 270 


5, 


fin 

.oU 


it 


0, 


,06 


7 


.92 + 0, 


,10 


1 


.U4 


n 


1 "7 


JUZio4 i , 


.oUH-UOZDlo. 


,7 


19, 


,954 + 


0, 


,027 


0, 


o 1 y 


it u. 


,069 


-0, 


,204 + 0, 


,040 


13510 ± 180 


7, 


no 

.uy 


it 


0, 


,08 


9 


.91 + 0, 


,18 


1 


no 

.UZ 


n 
U. 


Qn 


JU4zio4, 


CiA 1 CQnoc;i 


,7 


19, 


,789 ± 


0, 


,027 


0, 


fil c 


it U. 


,126 


-0, 


,216 + 0, 


,039 


11870 ± 200 


7, 


.y4 


1 

it 


0, 


,07 


11 


.57 + 0, 


,11 


U. 


A A 
.44 


U. 


,o ( 


JU ( UziD. 


oi 1111 nnn 


,0 


16, 


,061 ± 


0, 


,015 


0, 


'7 a. A 
. ( D4 


it U. 


,022 


-0, 


125 + 0, 


,019 


8800 ± 600 


6, 


nn 
,UU 


1 

it 


0, 


12 


9 


.08 + 0, 


,40 


U. 


Of; 
ZO 


n 
U. 


1 fi 

. iO 


JU / 40li. 


,0DH-iy4yZD. 


,5 


16, 


,259 ± 


0, 


,008 


0, 


,004 


it U. 


,018 


-0, 


,063 + 0, 


,020 


8190 ± 550 


5, 




1 

it 


0, 


12 


8 


.00 + 0, 


,20 


U. 


A 

40 


U. 


1 fi 
,iO 


JU / 4D10. 


.ooH-oyzzUo. 


1 


16, 


,589 ± 


0, 


,018 


0, 


DO / 


it u. 


,025 


-0, 


130 + 0, 


,026 


12130 ± 400 


5, 


yo 


it 


0, 


,12 


8 


.44 ± 0, 


,39 


n 


.4o 


n 
U. 


1 7 


jUozyu4, 


. 1 0+0 ( UOio. 


,4 


19, 


,242 ± 


0, 


,017 


0, 


. oyo 


it u. 


,048 


-0, 


,246 + 0, 


,026 


13430 ± 230 


7, 


.4Z 


it 


0, 


,08 


10 


.48 ± 0, 


,14 


n 
U. 


.0 1 


n 
U. 


.t5o 


JUo4oZ0. 


nO-LQVI f^E^I 

.uy+o 1 iODi. 


,7 


19, 


,890 + 


0, 


,034 


0, 


. ( UZ 


it u. 


,081 


-0, 


150 + 0, 


,044 


12950 ± 280 


7, 


. ( y 


it 


0, 


,08 


11 


.19 + 0, 


,17 


n 


.00 


n 
U. 


c;n 

.OU 


juyuuoz. 


fi/i _i_noQ/i 1 
,U4H-UZo4iv5. 


,8 


17, 


,965 ± 


0, 


,023 


0, 


OvSy 


it u. 


,033 


-0, 


130 + 0, 


,029 


8220 ± 330 


5, 


. / o 


it 


0, 


,07 


8 


.00 + 0, 


,20 


n 


OQ 

yo 


n 


1 fi 

.10 


Trim QOfi 


.Uo+o ( OoUo. 


,7 


18, 


,607 + 


0, 


,017 


0, 


ODU 


it U. 


,032 


-0, 


,235 + 0, 


,034 


12060 ± 250 


7, 


OI 

yi 


1 

it 


0, 


,08 


11 


.53 + 0, 


,14 


U. 


Ofi 

ZO 


n 
U. 


.00 


juyooOo. 


.D0H-4iUyz ( . 


,4 


19, 


,589 ± 


0, 


,022 


0, 


.004 


it u. 


,053 


-0, 


,248 + 0, 


,030 


13040 ± 490 


7, 


,00 


it 


0, 


,08 


10 


.91 + 0, 


,16 


n 
U. 


.04 


n 
U. 


A 

.40 


T1 1 1 OI 

J 1 i izlD. 


QO_l_1 1 1 V/l 
,oZH-lii ( 40. 


,0 


16, 


,235 ± 


0, 


,017 


0, 


oZU 


it u. 


,034 


-0, 


138 + 0, 


,026 


9400 ± 490 


5, 


CI 

.ol 


it 


0, 


12 


8 


.38 + 0, 


,27 


n 


.O / 


n 


1 fl 
. iO 


T1 1 Qni 7 


.40H-OO000U. 


1 


19, 


,413 ± 


0, 


,029 


0, 


'7fi1 


it U. 


,062 


-0, 


191 + 0, 


,045 


12030 ± 180 


7, 


1 1 


1 

it 


0, 


,08 


10 


.32 + 0, 


,12 


U. 


DO 


U. 


Qn 
oU 


T1 1 Q'ZOQ 


A A \ ^ OQ1 nc; 
,44H-lZolUu. 


,9 


19, 


,122 ± 


0, 


,031 


0, 


004 


it U. 


,059 


-0, 


160 + 0, 


,036 


13000 ± 390 


7, 


70 

. /y 


1 

it 


0, 


,08 


11 


.19 + 0, 


,18 


U. 


QO 

oy 


U. 


c;n 
,oU 


T1 1 /I QriQ 


CQ 1 Qfil Q /I Q 

.ooH-oDio4o, 


,8 


19, 


,678 ± 


0, 


,021 


0, 


OD ( 


it U. 


,057 


-0, 


,225 + 0, 


,035 


15090 ± 190 


7, 


oz 


1 

it 


0, 


,08 


10 


.54 + 0, 


,15 


U. 


D ( 


U. 


A fi 

.40 


T1 O'iQI a 


, iyH-iDUZU4, 


,6=^ 


19, 


,829 ± 


0, 


,018 


0, 


oUy 


it u. 


,068 


-0, 


,233 + 0, 


,028 


10920 ± 190 


5, 


1 o 

. iZ 


it 


0, 


,07 


8 


.00 + 0, 


,20 


z, 


QO 

.oZ 


n 
U. 


1 7 


T1 OQc;OQ 


. /oH-4/ouzy, 


1 


19, 


,291 ± 


0, 


,018 


0, 


iDOU 


it u. 


,080 


-0, 


136 + 0, 


,024 


10600 ± 200 


7, 


.Do 


it 


0, 


,08 


11, 


.45 ± 0, 


,13 


n 


.O / 


n 
U. 


Ad 
,40 


11 AA'\A9 

U Xl'iOrtZr. 


74-1-1 '^nQ'^8 


,6 


18, 


,578 ± 


0, 


,017 


0, 


757 


lb 0. 


,031 


-0, 


,060 + 0, 


,023 


8810 ± 220 


6, 


32 




0, 


,07 


9 


.67 + 0, 


,23 


0. 


61 


0. 


.17 


•J J- tj -LOiiU. 




,2 


17, 


,456 ± 


0, 


,019 


0, 


,687 


± 0, 


,026 


-0, 


,076 ± 0, 


,024 


9810 ± 320 


6, 


.66 




0, 


,06 


10 


.08 ± 0, 


,18 


0. 


.30 


0. 


.20 


11 ^21 22 




1 


18, 


,439 ± 


0, 


,017 


0, 


,794 


± 0, 


,041 


-0, 


190 + 0, 


,022 


15210 ± 250 


5, 


.60 




0, 


,08 


8 


.00 + 0, 


,20 


I 


.22 


0. 


.17 


11 "12(1^1 




,3 


18, 


,735 ± 


0, 


,018 


0, 


,789 


± 0, 


,041 


-0, 


,248 + 0, 


,028 


10340 ± 200 


5, 


.64 




0, 


,08 


8 


.00 + 0, 


,20 


I 


.40 


0. 


.17 


J153300, 


,03+492948. 


,3 


19, 


,221 ± 


0, 


,021 


0, 


,610 


± 0, 


,047 


-0, 


,247 + 0, 


,030 


15830 ± 320 


7, 


.75 




0, 


,08 


10 


.68 + 0, 


,20 


0, 


.51 


0, 


.49 


J161431, 


,28+191219, 


,4 


16, 


,187 + 


0, 


,019 


0, 


,893 


±0, 


,025 


-0, 


123 + 0, 


,024 


8590 ± 540 


5, 


.64 




0, 


12 


8 


.00 + 0, 


,20 


0, 


.43 


0. 


.16 


J161722, 


,51+131018, 


,8 


18, 


,605 ± 


0, 


,014 


0, 


,868 


±0, 


,041 


-0, 


189 + 0, 


,021 


10480 ± 240 


5, 


.78 




0, 


,08 


8 


.00 + 0, 


,20 


1, 


.32 


0. 


.17 


J174140, 


,49+652638, 


,7 


18, 


,271 ± 


0, 


,022 


0, 


,879 


±0, 


,047 


-0, 


177 + 0, 


,030 


9790 ± 240 


5, 


.19 




0, 


,06 


8 


.00 ± 0, 


,20 


1, 


.13 


0. 


.16 


J184037, 


,78+642312, 


,3 


18, 


,757 + 


0, 


,014 


0, 


,818 


±0, 


,047 


-0, 


108 + 0, 


,019 


9140 ± 170 


6, 


.16 




0, 


,06 


9 


.25 + 0, 


,13 


0, 


.80 


0, 


.17 


J191311, 


,59+372631, 


,7 


15, 


,430 + 


0, 


,007 


0, 


,722 


±0, 


,016 


+0, 


,030 + 0, 


,008 


8070 ± 630 


7, 


.04 




0, 


12 


11 


.55 + 0, 


,40 


0, 


.06 


0. 


.25 


J211921, 


,96-001825, 


,8=^ 


20, 


,000 + 


0, 


,021 


0, 


,867 


±0, 


,092 


-0, 


197 + 0, 


,033 


10360 ± 250 


5, 


.36 




0, 


,07 


8 


.00 + 0, 


,20 


2, 


.51 


0. 


.17 


J213513, 


,09-072442, 


,5 


19, 


,533 ± 


0, 


,020 


0, 


,584 


±0, 


,060 


-0, 


,212 + 0, 


,030 


11230 ± 260 


7, 


.30 




0, 


,08 


10 


.70 ± 0, 


,16 


0, 


.58 


0, 


.33 


J221928, 


,48+120418, 


,6 


17, 


,681 ± 


0, 


,017 


0, 


,915 


±0, 


,031 


-0, 


150 + 0, 


,022 


9700 ± 160 


5, 


.12 




0, 


,08 


8 


.00 + 0, 


,20 


0, 


.86 


0. 


.16 


J221936, 


,32-092617, 


,4 


19, 


,887 + 


0, 


,036 


0, 


,593 


±0, 


,126 


-0, 


,223 + 0, 


,053 


14490 ± 250 


7, 


.83 




0, 


,08 


11, 


.03 ± 0, 


,12 


0, 


.59 


0. 


.51 


J222859, 


,93+362359, 


,6 


16, 


,522 ± 


0, 


,011 


0, 


,674 


±0, 


,019 


-0, 


,025 + 0, 


,014 


8840 ± 540 


6, 


.73 




0, 


12 


10 


.63 + 0, 


,37 


0, 


.15 


0. 


.20 


J223630, 


,08+223223, 


,8" 


17, 


,004 + 


0, 


,021 


0, 


,663 


±0, 


,028 


-0, 


,211 + 0, 


,026 


11290 + 50 


6, 


.30 




0, 


,02 


8 


.83 + 0, 


,27 


0, 


.43 


0, 


.17 


J231757, 


,41+060252, 


1 


19, 


,449 ± 


0, 


,035 


0, 


,647 


±0, 


,065 


-0, 


104 + 0, 


,040 


12010 ± 280 


6, 


.98 




0, 


,08 


10 


.03 ± 0, 


,16 


0, 


.76 


0, 


.26 


J231910, 


,03+175824, 


,6 


19, 


,448 ± 


0, 


,017 


0, 


,610 


±0, 


,052 


-0, 


145 + 0, 


,024 


13840 ± 190 


7, 


.19 




0, 


,08 


10 


.02 + 0, 


,14 


0, 


.77 


0. 


.32 


J233705, 


,02+153353, 


,6 


19, 


,676 ± 


0, 


,022 


0, 


,614 


±0, 


,063 


-0, 


,224 + 0, 


,033 


11700 ± 180 


7, 


.63 




0, 


,07 


11 


.11 + 0, 


,11 


0, 


.52 


0, 


.45 



Brown ot al. (2010) 
^ .Kilic ct al.. t200a') 

using WD model spectra with atmospheric parameters 
customized for each target. The results are consistent 
within 10 km s~^, which we take as our systematic ve- 
locity uncertainty. The Appendix data table presents the 
full set of radial velocity measurements for the 7 newly 
identified ELM WDs with multiple observations. 

Six of the 7 newly identified ELM WDs with multiple 
observations display significant radial velocity variability. 
The ELM WD with no significant radial velocity variabil- 
ity is J0900-I-0234, however we cannot rule out whether 
it is a binary or not. We first observed J0900+0234 with 
a single exposure on 2011 March 3 and it had a heliocen- 
tric radial velocity of 64± 12 km s"^ On 2011 May 9 we 
re-observed the ELM WD with 9 back-to-back 2.5 min 
exposures and it had a mean 67 km s~^ velocity with a 
±24 km dispersion. Although the measurement error 
was ±13 km s~^, half of the observed dispersion, there 
was no obvious periodicity. Thus J0900+0234 shows no 
significant velocity variation over the observed time base- 
lines. A future epoch of observations is required to rule 
out the possibility that our existing observations may 
have sampled the same orbital phase. Follow-up obser- 



vations are planned for the other ELM WDs as well. 
2.6. Orbital Elements 

We now compute the orbital period and other orbital 
elements for the six ELM WDs with significant radial 
velocity variability. We begin by solving for the best-fit 
period that minimizes for ^ circular orbit. Figure HI 
plots the periodograms. In some cases there are multiple 
period solutions because of insufficient coverage, however 
in all cases the periods are constrained to be <1 day. We 
estimate the period error by conservatively identifying 
the range of periods with < "^Xmim where Xmin ^he 
minimum x^. 

We compute best- f it orb ital elements using the code 
of iKenvon fc Garcial (|1986[ ). which weights each velocity 
measurement by its associated error. The uncertainties 
in the orbital elements are derived from the covariance 
matrix and x^. To verify these uncertainty estimates, 
we perform a Monte Carlo analysis where we replace the 
measured radial velocity v with v + gSv, where 6v is the 
error in v and g is a Gaussian deviate with zero mean 
and unit variance. For each of 10000 sets of modified ra- 
dial velocities, we repeat the periodogram analysis and 
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TABLE 2 

Binary Orbital Parameters 



Object 


P 


K 


7 


Spec. Conjunction 


Mass Function 


Ah 


'J'mcrgc 




(days) 


(km s~^) 


(km s~^) 


(days + 2455250) 




(Mq) 


(Gyr) 


J011210.25+183503.7 


0.14698 ± 0.00003 


295 ± 2 


-121 + 1 


259.61940 + 0.00015 


0.392 + 0.007 


> 0.62 


< 2.7 


J015213. 77+074913.9 


0.32288 ± 0.00014 


217 + 2 


-61 + 1 


261.50244 + 0.00040 


0.341 + 0.008 


> 0.57 


< 22 


J090052. 04+023413.8^ 




< 24 


67 + 8 










J144342. 74+150938.6 


0.19053 ± 0.02402 


307 + 3 


-172 + 3 


490.77780 + 0.00057 


0.569 + 0.074 


> 0.83 


< 4.1 


J151826. 68+065813.2 


0.60935 ± 0.00004 


172 + 2 


-67 + 2 


25.98375 + 0.00131 


0.322 + 0.005 


> 0.58 


< 101 


J174140.49+652638.7 


0.06111 ± 0.00001 


508 + 4 


-70 + 3 


259.54392 + 0.00006 


0.830 + 0.018 


> 1.10 


< 0.17 


J184037.78+642312.3 


0.19130 ± 0.00005 


272 + 2 


-76 + 2 


259.43724 + 0.00086 


0.399 + 0.009 


> 0.64 


< 5.0 



Existing measurements are eonsistent witli no variation. 
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Fig. 4. — Periodograms for the 6 newly Identified ELM WDs with 
significant velocity variability. Some objects are well constrained, 
while others have multiple period aliases. In all cases the periods 
are <1 day. 

derive new orbital elements. We adopt the inter-quartile 
range in the period and orbital elements as the uncer- 
tainty. For binaries with multiple period aliases, both 
approaches yield similar uncertainties. When there are 
several equally plausible periods, the Monte Carlo anal- 
ysis selects all possible periods and derives very large 
uncertainties. In these cases, we adopt errors from the 
covariance matrix for the lowest orbital period. We 
plot the best- fit orbits in Figure [S] 

Table [2] presents the best-fit orbital parameters. 
Columns include orbital period (P) , radial velocity semi- 
amplitude (K), systemic velocity (7), the time of spec- 
troscopic conjunction (the time when the object is closest 
to us), mass function (see Eqn. 1 below), and minimum 
secondary mass (assuming i = 90°). The systemic ve- 
locities in Table [2] are not corrected for the WDs' grav- 
itational redshifts, which should be subtracted from the 
observed velocities to find the true systemic velocities. 
This correction is a few km s^^ for a 0.17 Mq helium 
WD, comparable to the systemic velocity uncertainty. 

3. RESULTS 

The orbital solutions constrain the mass and thus the 
nature of the ELM WD binary companions, as well as 
the binary systems' gravitational wave merger times. We 
discuss each binary in turn. 



3.1. J011210. 25+183503. 7 

The ELM WD J0112-f 1835 has a well-constrained or- 
bital period of 3.5275 ± 0.0006 hr and a radial velocity 
amplitude of 590 ± 4 km s^^. Its binary mass function is 
given by 



M| sin''i 
{Ml + A/2)' 



PK^ 
2ttG 



= 0.392 ± 0.0G7Mf 



0; 



(2) 



where i is the orbital inclination an gle. Mi c± 1 6 Mp, 
is the ELM WD mass inferred from lPanei et all (|2007[ i 
tracks, and M2 is the companion mass. For an edge-on 
orbit with i = 90°, Eqn. 2 provides the minimum com- 
panion mass (see Table ^ . Assuming a random orbital 
inclination distribution, on the other hand, allows us to 
calculate the probability of different companion masses. 

Given the observed orbital parameters, there is a 71% 
probability that J0112+1835's unseen companion is a 
WD with <1.4 Mq and a 14% probability that the com- 
panion is a neutron star with 1.4-3.0 Mq. The likelihood 
that the system contains a pair of WDs whose total mass 
exceeds the Chandrasekhar mass is 4%. If we assume 
the mean inclination angle for a random stellar sample, 
i = 60°, we get an estimate of the most probable compan- 
ion mass. For J0112-I-1835, the most likely companion is 
a 0.85 Mq WD at an orbital separation of 1.2 Rq. 

There is no evidence for a 0.85 Mq WD in the spectrum 
of JOl 124-1835, nor do we except there to be. If we pes- 
simistically assume that the two WDs formed at the same 
time 100 Myr - 1 Gyr ago, we would exp ect the 0.85 Mq 
comp anion to have Mg ^ 11 — 13 mag (jBergeron et al.l 
Il995f ): it would be 15 - 100 times less luminous than the 
0.16 Mq WD. Of course to form a short-period ELM WD 
binary like J0112+1835 requires two consecutive phases 
of common-envelope evolution in which the ELM WD is 
created last, giving the more massive secondary yet more 
time to cool and fade. 

To understand the evolutionary history of J0112+1835 
and our other ELM WD binaries requires assumptions 
about the energy balance and angular momentum bal- 
ance o f the common envel ope phase (e.g. lNelemans et alj 
[200I). IKilic et al.l (I2010b l) discuss one possible origin for 
a 1 h orbital period ELM WD evolving from a system 
containing a 3 Mq and a 1 Mq star. The 3 Mq star 
evolves off the main sequence, overfiows its Roche lobe 
as a giant with a 0.6 Mq core, forms a helium star (sdB) 
which does not expand after He-exhaustion in the core, 
and turns into a WD. The 1 Mq star also overflows its 
Roche lobe after main-sequence evolution when its core 
is around 0.2 Mq. We can estimate orbital separations if 
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we assume that the evolved stars exactly fill their Roche 
lobes. In this case, the first common-envelope phase has 
an orbital separation of 860 Rq and the second common- 
envelope phase has an orbital separation of 25 Rq. The 
orbital separation of J0112-I-1835 and our other ELM 
WD binaries is now ~1 Rq. 

General relativity predicts that short period binaries 
like J0112-I-1835 lose energy and angular momentum to 
gravitational wave radiation. The time scale for the bi- 
nary to shrink and begin mass transfer via Roche-lobe 
overflow is given by the gravitational wave merger time 

,^iMl±^pm^,0-^Gyr (3) 
Ml Ma 

where the masses are in M q and the period P is in hours 
(|Landau fc Lifshitd 119581 ). Inserting the minimum com- 
panion mass yields the maximum merger time given in 
Table [2] For the most probable companion mass of 0.85 
Mq, J0112+1835 will begin mass transfer in 2.1 Gyr. 
iKilic et all ()2010br ) discuss the many possible stellar evo- 
lution paths for such a system. This system's mass ratio 
M1/ M2 <0.26 suggest s that mass transfer will be sta- 
ble (jMarsh et al.ll2004[) and that J0112-K1835 wiU likely 
evolve into an AM CVn system. 

3.2. J015813. 77+074913. 9 

The ELM WD J0152-h0749 has a longer orbital period 
of 7.749 ± 0.003 hr and a radial velocity amplitude of 
434 ±4 km s"^ There is a 74% probability that the 
companion is a WD with <1.4 Mq and a 13% probability 
that the companion is a neutron star with 1.4-3.0 Mq. 
For i — 60°, the most likely companion is a 0.78 Mq WD 
at an orbital separation of 1.9 Rq. This system will not 
begin mass transfer within a Hubble time. 

3.3. JI44342.74+I50938.6 

The ELM WD J1443-H1509 has a best-fit orbital pe- 
riod of 4.573 hr. However, the current data set (which 
spans only 3 nights) allows for a significant alias at 5.75 
hr. The relatively large 614 ± 6 km s~^ radial veloc- 
ity amplitude of this system implies that the companion 
must be relatively massive, regardless of the exact pe- 
riod. Adopting the best-fit orbital period, there is a 60% 
probability that the companion is a WD with <1.4 Mq 
and a 20% probability that the companion is a neutron 
star with 1.4-3.0 Mq. For i = 60°, the most likely com- 
panion is a 1.15 Mq WD at an orbital separation of 1.5 

Rq. 

This system will begin mass transfer in less than 4.1 
Gyr. The likelihood that the system contains a pair of 
WDs whose total mass exceeds the Chandrasekhar mass 
is 6%. Given the observed mass ratio M1/M2 <0.20, this 
system will undergo stable mass transfer and will likely 
evolve into an AM CVn system. 

3.4. J151826. 68+065813. 2 

The ELM WD J1518-H0658 is similar to J0152-^0749 
except that the secondary is likely at a larger orbital 
separation. The system has an orbital period of 14.624 ± 
0.001 hr and a radial velocity amplitude of 344 ± 4 km 
s^^. Given these parameters, there is a 74% probability 
that the companion is a WD with <1.4 Mq and a 13% 
probability that the companion is a neutron star with 



1.4-3.0 Mq. For i = 60°, the most likely companion is 
a 0.78 Mq WD at an orbital separation of 3.0 Rq. This 
system will not begin mass transfer within a Hubble time. 

3.5. J174140.49+652638.7 

The ELM WD J1741-H6526 is arguably the most inter- 
esting of the six new systems. J1714+6526 has an orbital 
period of 1.4666 ± 0.0001 hr and a radial velocity ampli- 
tude of 986 ± 4 km s^^. Because our 6 minute exposure 
times span 7% of its orbital phase (^0 = 0.43 radians), 
the observed amplitude is underestimated by a factor of 
amS(j)/S(j) = 0.97. The true radial velocity amplitude of 
the ELM WD is thus 1,016 km s-\ 

Using the corrected orbital parameters, there is a 43% 
probability that J1741-|-6526's companion is a WD with 
<1.4 Mq and a 31% probability that the companion is 
a neutron star with 1.4-3.0 Mq. For i = 60°, the most 
likely companion is a 1.55 Mq neutron star at an orbital 
separation of 0.8 -Rq. Given that this putative neutron 
star must have accreted material from the common en- 
velope evolution of the ELM WD progenitor, a neutron 
star companion is possibly a milli-second pulsar. 

Milli-second pulsars in short-period orbits are diffi- 
cult to detect because of their rapidly changing velocity, 
but are valuable probes of general relativity and grav- 
itational wave physics. For the most likely companion 
mass of 1.55 Mq, J1741-h6526 wiU merge in 130 Myr 
- twice as fast as the Hulse- Taylor pulsar. The grav- 
itational wave strain for this system /i ~ 7 x 10~^^ is 
in principle detectable by the proposed LISA mission, 
but its orbital frequency places the system below the 
expected confusion-limit fr om other double-d egenerate 
gravitational wave sources (IRoelofs et al.ll2007l) . We are 
pursuing follow-up radio and X-ray observations to test 
whether or not J1741-I-6526 is a milli-second pulsar bi- 
nary. 

If the companion is a massive WD, on the other hand, 
there is a 30% likelihood that J1741-I-6526 contains a pair 
of WDs whose total mass exceeds the Chandrasekhar 
mass. With a M1/M2 <0.15 mass ratio, J1741-h6526 
will initially evolve into a stable mass-transfer AM CVn 
system. When the massive WD accretes sufficient mass, 
it is possible J1741-I-6526 will explode as a Type la su- 
pernova. If the companion is not massive enough to be 
a la progenitor, the s ystem will likely crea te an under lu- 
minous .la explosion (jBildsten et al l [20071 ). We are pur- 
suing time series photometry of J1741-I-6526 (Hermes et 
al., in prep.) to better constrain the orbital inclination 
and future evolution of this system. 

3.6. J184037.78+642312.3 

The ELM WD J1840-I-6423 has a 4.5912 ± 0.0012 hr 
orbital period, with a significant alias at 3.85 hr, and a 
544 ± 4 km s~^ radial velocity amplitude. Assuming the 
best-fit orbital period, there is a 70% probability that the 
companion is a WD with <1.4 Mq and a 15% probability 
that the companion is a neutron star with 1.4-3.0 Mq. 
For i — 60°, the most likely companion is a 0.88 Mq WD 
at an orbital separation of 1.4 Rq. 

This system is simflar to J1443-hl509: J1840-f 6423 wiU 
begin mass transfer in less than 5.0 Gyr, and the likeli- 
hood that the system contains a pair of WDs whose total 
mass exceeds the Chandrasekhar mass is 4%. Given the 
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Fig. 6. — Binary orbital period versus total system mass for the 
full ELM Sur vey and for the Supernova Progenitor Survey (SPY, 
IKoester et al . 2009). We plot total system mass assuming i = 
60° when orbital inclination is unknown, and the correct system 
mass when inclination is known. WD binaries from SPY are drawn 
with open triangles, previously published ELM Survey binaries are 
drawn with solid red circles, and the 6 new ELM WD binaries from 
this paper are drawn with solid blue diamonds. The dashed line 
shows the threshold at which a 1:1 mass ratio system will merge in 
10 Gyr. 

observed mass ratio M1/M2 <0.27, this system should 
undergo stable mass transfer and will likely evolve into 
an AM CVn system. 

4. DISCUSSION 

The ELM Survey has now identified 19 merging ELM 
WD systems that will coalesce in <10 G yr. The first 
12 were summarized bv lKilic et al.l (|2011aD . Three more 
systems were published earlier this year: two 39 minute 
orbital period binaries (Kilic et al. 2011c. b) and one 12 
minute period eclipsing binary (Bro wn et al]i20ilb[ ). In 
this paper we add four more merging systems to the 
count. To aid the reader, Table [3] summarizes the prop- 
erties of the 19 merging ELM WD systems. 

While the merging ELM WD sample is not complete, 
properties such as orbital period and secondary mass 
are in principle independent of color and magnitude. 
Thus we split the merging ELM WD sample into thirds 
and search for significant correlations among the sample 
properties. 

Looking at Tabled the six ELM WDs with <1.1 hr 
orbital periods are 6000±2400 K hotter than the seven 
ELM WDs with >3.5 hr orbital periods. In other words, 
we observe an absence of cool ELM WDs with short 
orbital periods. This period-temperature dependence 
makes sense if the ELM WDs in short orbital period sys- 
tems merge before they cool0 We take this correlation as 
evidence that common envelope evolution creates ELM 
WDs directly in short orbital period systems. 

The minimum companion masses of the ELM WDs 
span an order of magnitude, from 0.1 Mq to 1.1 Mq, and 
also appear to correlate with ELM WDs' temperatures. 

^ The period- temperature elfect is unchanged if we drop the 12 
minute orbital period syste m, which is pos sibly heating up due to 
tidal effects tPiro 2011: F uller fc Lai|[20TlD . 



The coolest seven ELM WDs, those with roff< 10,500 K, 
have minimum companion masses that are 0.41±0.23 Mq 
larger than the hottest six ELM WDs with roff>16,000 
K. This result largely comes from our present survey, 
which targets relatively cool ELM WDs and finds only 
extreme mass-ratio binaries. Our interpretation is that 
shorter period binary systems experience increased mass 
loss during their evolution and end up with lower mass 
companions. A larger sample is required to establish 
these trends with increased significance. 

Figure [S] compares the distribution of binary orbital 
period and total system mass for the ELM Survey 
with that of the Supernova Progenitor Su rvey (SPY, 
iNapiwotzki eral|[200lt [Koester et all [20091 ) . The SPY 
survey has measured the velocity variability of ~1000 
WDs, the lar gest survey o f its kind. Our Figure is an 
adaptation of 'Geier etall (|2010l )'s Fi gure 4, where we 
plot total system mass assuming i = 60° when orbital in- 
clination is unknown, and the correct system mass when 
inclination is known via ellipsoidal variations and/or 
eclipses. It is notable that the SPY survey, which sam- 
ples the full WD population, finds only a handful of sys- 
tems with orbital periods short enough that they might 
possibly merge in less than 10 Gyr. 

The ELM survey, which targets '-0.2 Mq WDs, finds 
almost every object in a system with <1 day orbital pe- 
riod, the majority of which will merge due to gravita- 
tional wave radiation in less than 10 Gyr. The merg- 
ing ELM WD systems are unlikely Type la supernovae 
progenitors, however, because the total mass of the sys- 
tems is likely below the Chandrasekhar mass. Their most 
likely evolutionary futures include the formation of sta- 
ble mass-transfer AM CVn systems and undcrluminous 
supernovae, or unstable mass-transfer mergers that form 
extreme helium stars (RCrB) and sin gle helium-enriche d 
subdwarf O stars (discussed further in lKilic et al.l2010b| l . 
One approach to constrain these scenarios is to com- 
pare ELM WD merger rates with the form ation rates of 
different classes of objects as attempted in iBrown et al.l 
([201 la). A larger, well-defined sample will provide im- 
proved constraints, as will follow-up light curves that di- 
rectly measure the orbital inclination and nature of the 
unseen binary companions. 

We close by noting that at least one-third of the merg- 
ing ELM WD systems have the kinematics and locations 
of halo objects. Looking at Tabled four ELM WD sys- 
tems (J0112, J0818, J1443, NLTT 11748) have systemic 
radial velocities I7I > 100 km s^^. Proper motions with 
5 mas y r~^ uncertainties are available for 15 of the 19 
objects (jMunn et al.ll2004HLepine fc Sharall2005l ). Com- 
bining radial velocities and proper motions reveals five 
systems (J0106, J0818, J1053, J1443, NLTT 11748) with 
total space velocities >200 km s~^ with respect to the 
Sun. Thus 6 unique ELM WD systems (32%) have kine- 
matics that indicate a halo origin. In addition to mo- 
tions, physical locations also support a halo origin. The 
ELM WDs are not clustered at low Galactic latitudes in 
our survey, as expected for a disk population, but rather 
are found equally at high and low Galactic latitudes. The 
typical ELM WD in our survey has median luminosity 
Mg = 8.5, apparent magnitude 50 = 18.8, and Galactic 
latitude b = 43°, and thus is located ~1 kpc above the 
Galactic plane. We conclude that ELM WD systems con- 
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TABLE 3 

Merger Systems in the ELM Survey 



Object 




log 9 


Mass 


P 


K 


M2 


''"merge 


7 




MDcc 


Rcf 




(K) 




Mq 


days 


km s~^ 


Mq 


Gyr 


km s~^ 


mas yr~^ 


mas yr~^ 




JUUzz— 11)14 


18980 


7.15 


0.33 


0. 07989 


145.6 


> 0.19 


< 0.73 


81 


—7.8 


— 13.2 


5 


JUIUd— lUUU 


16490 


6.01 


0.17 


0.U2715 


395.2 


0.43 


0.037 


2 


20.2 


—9.6 


6 


JUli2+lo35 


9770 


5.57 


0.16 


0.14698 


295.3 


> 0.62 


< 2.67 


— 121 


7.8 


— 17.4 





JUd51+2o44 


16400 


6.79 


0.25 


0.00885 


657.3 


0.55 


r\ r\r\r\r\ 

U.UUU9 


17 


—3.6 


— 1.2 


2 


JU755+4y(jD 


13160 


5.84 


0.17 


U.UboUz 


438.0 


> 0.81 


< 0.22 


—51 






1 


JUolo+oooo 


10620 


5.69 


0.17 


0.18315 


170.0 


> 0.26 


< 8.89 


— 201 






1 




8880 


6.44 


0.17 


r\ o /I ^ r\r\ 

U.Z44UU 


271.1 


> 0.76 


< 8.42 


—52 


3.4 


— 19.2 


3 




1 0900 


u.zo 


17 


n 07870 


ouu.y 




< n 47 




7 


A 


Q 




J0923+3028 


18350 


6.63 


0.23 


0.04495 


296.0 


> 0.34 


< 0.13 


2 


-4.2 


-24.9 


1 


J1053+5200 


15180 


6.55 


0.20 


0.04256 


264.0 


> 0.26 


< 0.16 


12 


-29.7 


-31.2 


3,8 


J1233+1602 


10920 


5.12 


0.17 


0.15090 


336.0 


> 0.86 


< 2.14 


-35 






1 


J1234-0228 


18000 


6.64 


0.23 


0.09143 


94.0 


> 0.09 


< 2.69 


94 


-14.5 


-12.3 


5 


J1436+5010 


16550 


6.69 


0.24 


0.04580 


347.4 


> 0.46 


< 0.10 


-30 


7.8 


-5.1 


3,8 


J1443+1509 


14770 


6.06 


0.17 


0.19053 


306.7 


> 0.83 


< 4.09 


-172 


-30.9 


-53.5 





J1630+4233 


14670 


7.05 


0.30 


0.02766 


295.9 


> 0.30 


< 0.03 


-11 


2.3 


-7.3 


7 


J1741+6526 


9900 


5.20 


0.16 


0.06111 


508.0 


> 1.09 


< 0.17 


-70 


-3.2 


-4.6 





J1840+6423 


9100 


6.22 


0.17 


0.19130 


272.0 


> 0.64 


< 5.00 


-76 


-12.4 


-26.7 





J2119-0018 


10360 


5.36 


0.17 


0.08677 


383.0 


> 0.75 


< 0.54 


-28 






1 


NLTT 11748 


8690 


6.54 


0.18 


0.23503 


273.4 


0.76 


7.20 


126 


236.1 


-179.2 


4,9,10 



Referenc es. — (0) this pa per; (1) IBrown et all II2010I) : (2) [ Brown et aU i2011t): (3') IKilic et al.l ll2010bl) : (4) IKilic et al 
Kilic et al.l Il2011ari : (e VlKi ic et all I l201lj) r77') iKiUc et al.l I l2011bl) : f8) (Mullallv et all 1 120091) ISteinfadt et al 



2010al ): (5) 
1 20101) : (10) 



Kawka et al.l II2010I) 



Note. — Measurement errors reported in the references. Proper motions are from USNO-B+SDSS ||M unn et al.|[200l) . 



tinue to form and merge in both the disk and the halo, 
a conclusion that h as implications for gravitational wave 
source predictions (jRuiter et al.ll2009l ). 



5. CONCLUSION 

We present a targeted spectroscopic survey of ELM 
WDs candidates selected by color. The survey is suc- 
cessful: it is now 71% complete and has uncovered 18 
new ELM WDs. Of the 7 ELM WDs with follow-up 
observations, 6 are compact binary systems and 4 have 
gravitational wave merger times less than 5 Gyr. The 
most intriguing new object is J1741-)-6526, which likely 
has either a milli-second pulsar binary companion or a 
massive WD companion making the system a possible 
supernova Type la or .la progenitor. Follow-up observa- 
tions are underway to establish the nature of this system 
as well as the other ELM WDs. 

Based on these initial results, we expect that com- 
pleting our targeted ELM WD survey will double the 
number of merging systems in our sample from 19 to 
~40 systems. We expect that photometric follow-up will 
reveal additional eclipsing systems. Our efficiency for 
ELM WDs discoveries increases with apparent magni- 
tude such that, if we were to expand our discovery sur- 
vey to g ~ 20.5, we could in principle double again our 



sample of ELM WDs (although the observations are more 
expensive). The absence of short- merger time systems at 
cooler temperatures suggests that expanding our survey 
to hotter objects may yield additional ~10 min orbital 
period (<1 Myr merger time) systems. These are direc- 
tions we will pursue in upcoming observing runs. With 
a sample of ^100 ELM WDs systems spanning a well- 
defined range of temperature, we look forward to placing 
robust constraints on the role of these detached double 
degenerate binaries as supernovae progenitors and grav- 
itational wave sources. 

We thank M. Alegria, J. McAfee, and A. Milone for 
their assistance with observations obtained at the MMT 
Observatory, and P. Berlind and M. Calkins for their as- 
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is managed by the Astrophysical Research Consortium 
for the Participating Institutions. This research makes 
use of NASA's Astrophysics Data System Bibliographic 
Services. This work was supported in part by the Smith- 
sonian Institution. MK was supported in part by NASA 
through the Spitzer Space Telescope Fellowship Program, 
under an award from CalTech. 
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APPENDIX 
DATA TABLE 

Table 4 presents our radial velocity measurements. The Table columns include object name, heliocentric Julian 
date (based on UTC), heliocentric radial velocity (uncorrected for the WD gravitational redshift), and velocity error. 
[See on-line journal or tab4.dat in the arXiv submission.] 
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